Time domain Einstein-Podolsky-Rosen correlation 



Nobuyuki Takei 1,2 |j Noriyuki Lee 1 ' 2 , Daiki Moriyama 1 , J. S. Neergaard-Nielsen 1 ' 3 , and Akira Furusawa 1,2 
1 Department of Applied Physics, School of Engineering, The University of Tokyo, 
7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan 
2 CREST, Japan Science and Technology Agency, 1-9-9 Yaesu, Chuo-ku, Tokyo 103-0028, Japan 
3 Niels Bohr Institute, University of Copenhagen, DK 2100, Denmark 
(Dated: February 1, 2008) 

We experimentally demonstrate creation and characterization of Einstein-Podolsky-Rosen (EPR) 
correlation between optical beams in the time domain. The correlated beams are created with 
two independent continuous-wave optical parametric oscillators and a half beam splitter. We define 
temporal modes using a square temporal filter with duration T and make time-resolved measurement 
on the generated state. We observe the correlations between the relevant conjugate variables in time 
domain which correspond to the EPR correlation. Our scheme is extendable to continuous variable 
quantum teleportation of a non-Gaussian state defined in the time domain such as a Schrodinger 
cat-like state. 
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Quantum entanglement is one of the most fundamen- 
tal properties of quantum mechanics and also the essen- 
tial ingredient in quantum information processing 0, Q . 
The nonclassical correlation was first argued by Einstein, 
Podolsky, and Rosen (EPR) relating to canonically con- 
jugate continuous variables (CVs) Qj. This originally de- 
vised EPR correlation has been experimentally realized 
by employing a two-mode squeezed vacuum state 0, |^ , 
where the relevant variables are quadrature-phase am- 
plitudes of optical fields. Since this significant mile- 
stone, there have been further experimental observations 
of the EPR correlations 0, H @ an d theoretical char- 
acterizations ^JEl- Moreover the EPR correlation, or 
the two-mode squeezed vacuum has been used as a re- 
source in implementation of quantum protocols. Thanks 
to well-established tools of manipulating the EPR cor- 
relation, CV quantum information processing (QIP) has 
been rapidly developing in recent years and has been 
successful in implemen ting qu antum protocols such as 
quantum teleportation p^lla.ll4j . 

Quantum teleportation is one of the fundamental 
quantum protocols which exploit the EPR correlation 
or quantum entanglement 0, 0. By using an opti- 
cal field, it has been experimentally realized for Gaus- 
sian input states; a coherent state 0,0], a squeezed 
state Jl6l - and an EPR state, i.e., entanglement swap- 
ping |14j. These successes are based on well-developed 
techniques of optical Gaussian operations consisting of 
beam splitters, phase shifting, squeezing, phase-space 
displacement and homodyne detection. In addition to 
the teleportation, CV quantum protocols have so far- 
been performed only with Gaussian states and opera- 
tions. However non-Gaussian states or non-Gaussian op- 
erations are further required to extract versatile potential 
of CV QIP 0]. Thus quantum teleportation of a non- 
Gaussian state would become the next important chal- 
lenge toward universal CV QIP. 

For the generation of a non-Gaussian state in an op- 



tical system, there is considerable interest in the tech- 
niques where photon counting is performed for a tapped- 
off beam from a Gaussian state light beam. When count- 
ing a photon, the tapped Gaussian state is converted 
to a non-Gaussian state, a technique called "photon- 
subtraction" . Recently there are some reports which 
show the generation of non-Gaussian states (Schrodinger 
cat-like states) with this technique fl8t |pf^ |. Since the 
photon counting event is defined in the time domain, con- 
ventional CV teleportation schemes with a sideband of an 
optical carrier |l3| cannot be used for the teleportation 
of such states. Therefore, in order to realize CV telepor- 
tation of such non-Gaussian states, one need to generate 
EPR correlation in the time domain. 

In this Letter we demonstrate experimental genera- 
tion and characterization of a two-mode squeezed vacuum 
state in the time domain, a step towards CV quantum 
teleportation of non-Gaussian states like a Schrodinger 
cat-like state. The main feature of our EPR resource is 
that we use two squeezed vacuum states of continuous- 
wave (CW) light beams from two independent subthresh- 
old optical parametric oscillators (OPOs) to obtain well- 
defined frequency and spatial modes, and also that we 
utilize almost the whole frequency bandwidth of the 
OPO cavities to define a quantum state in the time do- 
main. The use of the broad bandwidth allows us to make 
time-resolved measurements like photon counting as de- 
scribed below. The experiment presented here is quite 
different from the previous works in the frequency do- 
main 0, E| 0, which have dealt with only frequency 
sidebands a few MHz apart from the carrier frequency, 
and therefore the previous works are not compatible with 
photon counting. Our scheme would have advantages 
even over the time-domain pulsed scheme 9] with respect 
to spatial modes. The well-defined frequency and spatial 
modes are important for efficient scaling of a quantum 
circuit via interference of several beams. 

Moreover our scheme is compatible with CV entan- 
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glement distillation with the photon subtraction opera- 
tion j2jj. Combining the entanglement distillation and 
teleportation would result in the improvement of tele- 
portation fidelity |2l|. Our entangled state could also 
be applicable to loophole-free tests of Bell's inequali- 
ties [H |23|]. 

For an ideal two-mode squeezed vacuum state with in- 
finite squeezing, the EPR correlations can be described 
by the Wigner function of the two modes A and B: 

Wepr(za,Pa;:cb,Pb) = cS(x A - x b )S(pa +Pb), (!) 

where x,p are the quadrature amplitudes corresponding 
to the conjugate quadrature operators x,p which obey 
the commutation relation [x,p] = i/2 (H = 1/2) in anal- 
ogy with the position and momentum operators. This 
Wigner function expresses that the quadrature ampli- 
tudes of the two modes are perfectly correlated (xa = 
Xb) and anti-correlated (pa = ~Pb), respectively. In the 
real experiment where only finite squeezing is available, 
the state is no longer the ideal EPR state and shows finite 
variances in ([A(xa — £b)} 2 ) and ([A(pA + Pb)] 2 )- 

The two-mode squeezed vacuum generated by sub- 
threshold OPOs has EPR correlations within the vari- 
ance spectra S x (&), S P ({1) as usually measured in the 
frequency domain experiments. Here, f2 is a sideband 
frequency around the fundamental wavelength. When 
measuring the CW beams in the time domain we need to 
specify a certain temporal mode within which the entan- 
gled quantum state is defined. Because of the broadband 
correlations of the CW beams, there is a large amount of 
freedom in choosing a mode which will suit an intended 
application and will show significant correlation below 
the vacuum level. For instance, in a teleportation exper- 
iment the temporal mode could be chosen to match that 
of the input state. 

A simple choice of temporal mode is the square filter 
of duration T . Defining the filtered quadrature operators 
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and similarly for p\ these operators obey the same com- 
mutation relation, \x* ,p'\ — i/2. The EPR variances 
within this temporal mode will be given by |24j 
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<[A(p A +pb)?) = 5- / s p (n)-^dn. (4) 

The variances are thus filtered by a sine function, and 
by adjusting the integration time T the frequency range 




FIG. 1: Schematic setup of the experiment. OPOs: sub- 
threshold optical parametric oscillators, HBS: a half beam- 
splitter, PZTs: piezo-electric transducers, LOs: local oscilla- 
tors, ADCs: an analogue-to-digital converter. 



which contribute to the integrated correlation can be se- 
lected. The square temporal filter is simple, but not nec- 
essarily the most appropriate for an application of the 
time domain entanglement. Theoretical investigations 
concerning the shape of the filter combined with photon 
counting experiments have been done in Refs. |2J,|25j. 

The schematic diagram of our experiment is illustrated 
in Fig. 1. This is almost the same as the entanglement 
preparation part in our experiment of quantum telepor- 
tation The primary source of the experiment is a 
CW ThSapphire laser at 860nm, most of whose output 
is frequency doubled in an external cavity. The output 
beam at 430nm is divided into two beams to pump two 
OPOs. 

A two-mode squeezed vacuum is produced by combin- 
ing two squeezed vacua at a half beam splitter (HBS). 
Each squeezed vacuum is generated from a subthreshold 
OPO with a 10-mm-long KNb03 crystal. The crystal is 
temperature-tuned for type-I noncritical phase matching. 
Each OPO cavity is a bow-tie-type ring cavity consisting 
of two spherical mirrors (radius of curvature 50mm) and 
two flat mirrors. The round trip length is about 500mm 
and the waist size in the crystal is 20/im. An output 
coupler has transmissivity of 12.7%, while the other mir- 
rors are highly reflective coated at 860nm. They have 
a high transmission for 430nm so that the pump beam 
passes the crystal only once. The pump power is about 
70mW for each OPO. The total intracavity losses are 
around 2%, giving a cavity bandwidth of 7MHz HWHM. 
The resonant frequency of the OPO is locked via the 
FM sideband locking method by introducing a lock 
beam which is propagating against the squeezed vacuum 
beam to avoid the interference between the two beams. 
However a small fraction of the lock beam reflecting from 
surfaces of the crystal circulates backward and contam- 
inates the squeezed state. This problem is resolved by 
changing the transverse mode and the frequency of the 
lock beam, as described in detail in Ref. 



3 



DO 

;o 

<n 
CD 
o 
c 

"i_ 

> 




: (b) 









4 6 
Frequency (MHz) 



10 



FIG. 2: Calculated Fourier analysis of the 50M-sampled raw 
data without average, (a) S X (Q) for {[A(x A — £g)] 2 ). (b) 
S P (Q) for ([A(p£ +j5g)] 2 ). Each trace is normalized to the 
corresponding vacuum level. 




FIG. 3: Typical measured correlations. Only 50 points are 
picked up. (a) and (b) are measured quadrature values within 
time interval T = 0.2/isec for x and p quadratures, respec- 
tively. In each figure, trace (i) is for Alice, while (ii) for Bob. 



The output beams from the HBS are sent to Alice and 
Bob, and then measured using homodyne detectors with 
a bandwidth of 8.4 MHz. We lock relative phases between 
EPR beams and local oscillators (LOs) in the detection, 
i.e., x or p quadrature, in the following way. First, weak 
coherent beams which pass through the same paths as 
squeezed states are injected into the OPOs from one of 
the flat mirrors and used for the conventional dither and 
lock method. The relative phases between the weak co- 
herent beams and the squeezed vacua and between the 
weak coherent beams at the half beam splitter (HBS) 
are actively controlled by applying feedback voltages to 
piezos (PZTs). The quadratures to be detected in the 
homodyne detectors are also determined by locking the 
phases between the LOs and the injected coherent beams. 
However these weak beams and the modulation on them 
contaminate EPR correlation in some frequency ranges. 
Therefore we need to remove such beams and lock the 
phases without the beams. To do so, we introduce an 
electronic circuit that holds the feedback voltage and 
keeps the phase relation for a while. Within 2msec, the 
beams are blocked with mechanical shutters before the 
OPOs, and the EPR beams are measured with the ho- 
modyne detectors. Each output from the detectors is 
sampled with an analogue-to-digital converter (ADC) at 
a rate of 50M samples per second and then filtered on a 
PC to yield a number of measured quadrature values. 

Before observing the EPR correlation in the time do- 
main, we estimate what range of the frequency band- 
width can be used to yield a time-resolved quadrature 
value from the OPOs. Figure 2 shows the frequency 
spectra S x (£i) and S p (Cl) for the measured EPR beams 
calculated by performing Fourier transformation digitally 



on the 50M-sampled raw data. The EPR correlation is 
observed over the full bandwidth of the OPOs. fn par- 
ticular it is observed for low frequencies down to 5kHz (a 
high pass filter with a cut-off of 5kHz is used to eliminate 
noise at frequencies close to DC). The correlation at low 
frequencies is essential for the photon subtraction exper- 
iment for non-Gaussian operation or entanglement distil- 
lation |2sj . From the results, we can define a quantum 
state within a time interval of the inverse of the cavity 
bandwidth 7MHz. But the EPR correlations degrade at 
higher frequencies, so we use the temporal filter (2) with 
integration time T — 0.2/^sec, yielding 10000 points in 
each measurement. It follows from the frequency filters 
of sine fuctions in Eqs. (3) and (4) that we select mainly 
a frequency range below 5MHz. 

First we compare explicitly the time-resolved quadra- 
ture values measured at Alice and Bob in every 0.2^sec, 
rather than the variances of ([A(x^— ig)] 2 ) and ({A(p^ + 
Pb)] 2 ) performed in the previous works SJlIlli This 
is also different from the pulsed scheme |9( where EPR 
beams were recombincd at a beam splitter and returned 
to two squeezed vacua and then one of which was mea- 
sured in homodyne detection. Figure 3 shows the typ- 
ical example of the measured quadrature values within 
time interval T = 0.2^tsec, where only 50 points are 
picked up. As mentioned above, the measured values 
behave in such a way that x and p quadratures are corre- 
lated (x A — £ B ) and anticorrelated {jp A ~ — Pjj), respec- 
tively. This manifests the originally devised correlations 
discussed by EPR. 

We then estimate the degree of EPR correlation us- 
ing the correlation diagrams with 10000 points as shown 
in Fig. 4. These results show clear correlations in 



4 



(a) (b) 




-4 -2 2 4 -4 -2 2 4 

A A FA 



FIG. 4: Measured correlation diagrams of 10000 measured 
values for (a) x and (b) p quadratures. 

the relevant quadratures. By repeating the same mea- 
surement ten times, we obtained the correlations of 
{[A{x{-x f B )} 2 } = -3.30±0.28dB and ([A(p{+p f B )} 2 } = 
—3.74 ± 0.32dB compared to the corresponding vacuum 
variances. Accordingly the sufficient criteria [Hj, [llj for 
quantum entanglement holds; ([A(ar^ — x B )] 2 ) + ([A(p^ + 
Pb)} 2 ) = 0-45 ± 0.02 < 1. Therefore our generated 
state is entangled in the temporal mode defined within 
T = 0.2/isec. Note that we can decrease the integration 
time T for fast data acquisition and QIP, while the EPR 
correlation will degrade due to the contribution of higher 
frequency ranges. One would be able to obtain correla- 
tion over a broader band by use of small-size OPOs with 
shorter round-trip length or waveguide crystals, e.g., pe- 
riodically poled Lithium Niobate waveguides [29j ■ 

In summary we have experimentally demonstrated 
generation and characterization of a two-mode squeezed 
vacuum state in the time domain. We have observed 
the nonclassical correlation between the temporally mea- 
sured quadrature values. Our setup is almost the same as 
an entanglement preparation part of our quantum tele- 
portation setup With slight modification of classi- 
cal channels, we will be able to perform the broadband 
teleportation whereby it is possible to transfer quantum 
states defined in the time domain like a Schrodinger cat- 
like state, a single photon state and other non-Gaussian 
states. Moreover our scheme is compatible with CV en- 
tanglement distillation with photon subtraction [2fij ] . Our 
system has the well-defined frequency, spatial and tempo- 
ral modes. Therefore we could efficiently scale up quan- 
tum circuits by interfering several beams toward the uni- 
versal CV QIP. 

This work was partly supported by the MEXT and the 
MPHPT of Japan, and Research Foundation for Opto- 
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